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Abstract 

For the cylindrically symmetric ’’asymptotically flat” Einstein equations in the case 
of electro-vacuum it is known that solutions exist globally and also that this class of 
spacetimes is causally geodesically complete. Hence strong cosmic censorship holds for 
this class. An interesting question is whether these results can be generalized to include 
spacetimes with phenomenological matter, e.g. collisionless matter described by the 
Vlasov equation. Spherically symmetric asymptotically flat solutions of the Einstein- 
Vlasov system with small initial data are known to be causally geodesically complete. 
For arbitrary (in size) data it has been shown that if a singularity occurs, the first one 
occurs at the center of symmetry. In this paper we begin to study the question of global 
existence for the cylindrically symmetric Einstein-Vlasov system with general (in size) 
data and we show that if a singularity occurs at all, the first one occurs at the axis of 
symmetry. 


1 Introduction 

Consider a large ensemble of particles, e.g. stars in a galaxy, which interact by the 
gravitational field which they create collectively. In particular, assume that collisions 
are sufficiently rare to be neglected. In the Newtonian case one arrives at the so called 
Vlasov-Poisson system. From a mathematical point of view this system is quite well 
understood. The question of global existence for the Vlasov-Poisson system in three 
space dimensions for general smooth initial data was first solved by Pfaffelmoser HU. 
The relativistic analogue is much more complicated due to the fact that the field is 
governed by the Einstein equations instead of the Poisson equation. This is the so 
called Einstein-Vlasov system. For a detailed review concerning global results for the 
Einstein-Vlasov system, see [ 2 ] 

As is well known, freely falling particles follow the geodesics of spacetime. By assum¬ 
ing that collisions are sufficiently rare to be neglected, the matter distribution function 
/, i.e. a function / defined on phase space which describes the density of particles at po¬ 
sition x with momentum p at time t, is conserved along the geodesics. Hence / satisfies 
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the Vlasov equation 


df_ P v df 
dt p° dx v 


r u bc pV 


df_ 

dpV 


0, b,c = 0,1,2,3, v = 1,2,3, 


( 1 ) 


where are the Christoffel symbols associated to the metric in the local coordinates 
and p° is a solution to g a bP a P b = — 1- For more details concerning relativistic kinetic 
theory, see e.g. fE! or J], Now the Vlasov equation is coupled to the Einstein equations 


iy blfjab SvrT'jjfi, 


( 2 ) 


where the energy-momentum tensor T a b is given by 


Tab = -[ Iffl 2 dp. (3) 

Jr s Po 

Here t = x°, x = (x 1 ,^ 2 ,^ 3 ), p = (p 1 ,p 2 ,p 3 ), \g\ = det (g a b) , a,b = 0,1,2,3 and 
c = G = 1 where c is the velocity of light and G is the gravitational constant. 

When an isolated body is studied, the spacetime is assumed to be flat “far away” 
from the body. These spacetimes are called asymptotically flat. The topology of the 
initial surface is M 3 and the metric has to satisfy certain fall-off conditions. Intuitively 
the metric should tend to the Minkowski metric far away from the body in some coordi¬ 
nate system. This class is compatible only with a few different symmetry classes where 
spherical- and axial symmetry are the most important ones. In the class of asymptot¬ 
ically flat spacetimes only the spherically symmetric case has been considered. In m, 
na a global existence result for the spherically symmetric case with small initial data 
is obtained and it is proved that these spacetimes are geodesically complete. In PH it 
is shown, for arbitrary (in size) initial data, that if a singularity occurs, the first one 
occurs at the center of symmetry. The purpose of this paper is to prove a similar result 
for cylindrically symmetric spacetimes. 

Cylindrically symmetric spacetimes, admit asymptotic flatness conditions in two di¬ 
rections of the initial surface and is therefore not asymptotically flat. The reason that 
cylindrical spacetimes do not belong to the class of asymptotically flat spacetimes is 
that the symmetry group contains translation symmetry. However following Berger, 
Chrusciel and Moncrief [Hj we can state an asymptotic flatness condition in two of the 
three space dimensions on the initial surface (see section 2 below for a definition). For 
cylindrically symmetric, asymptotically flat electro-vacuum spacetimes it is known that 
global existence holds |H|, and in addition that these spacetimes are causally geodesically 
complete. This is done by rewriting the Einstein equations as a rotationally symmet¬ 
ric wavemap from M 1+2 into a Riemannian manifold. The target manifold satisfies a 
certain condition on the behaviour of the curvature at infinity and a convexity con¬ 
dition. Hence the results of Christodoulou and Tahvildar-Zadeh 0 0 apply. When 
matter is introduced this wavemap structure is lost. However it is possible to rewrite 
the Einstein equations as a wavemap with a forcing term from a rotationally symmetric 
2 + 1-dimensional curved spacetime with the same target manifold as in the vacuum 
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case. The forcing term corresponds to the matter term. Unfortunately we are not able 
to use this result and generalize the Christodoulou/Tahvildar-Zadeh result to this type 
of wavemap. Instead we use the method of Andreasson T], for cosmological spacetimes. 
In [lj there is no center of symmetry and by using that method we are not able to control 
the field or the matter at the axis of symmetry. Hence we need to introduce conditions 
both on the held and on the matter distribution function in a vicinity of the axis of 
symmetry to obtain global existence. However in the so called polarized case, i.e. when 
the metric is diagonal, we are able to remove one of the assumptions and still obtain 
the same result by keeping the condition on the matter. This is in analogy with the 
spherically symmetric case in m- The outline of the paper is as follows. In section 2 we 
define cylindrical spacetimes and the corresponding definition of asymptotic flatness for 
this case. In section 3 the main theorems are stated and the proofs of the theorems are 
outlined in section 4 and 5 respectively. Theorem 2 is proved by the strategy outlined 
in [Tj[. However to bound the matter terms we need to introduce a rescaled momentum. 
Theorem 1 is proved by Sobolev methods combined with the proof of Theorem 2. 

2 Cylindrical symmetry and asymptotic flatness 

By cylindrically symmetric initial data we mean that there exists a 2-dimensional isom¬ 
etry group G = U{ 1) x R, acting on a spacelike initial surface S = K 3 . The metric h 
and the second fundamental form K^ v of S satisfy 

£ x a h[iu — £x°-K^u — 0 ? 

where X a , a = 2,3 are two Killing vectors that generate the isometry group. Further¬ 
more the matter distribution function /o is also invariant under the action of G, i.e. 
fo(Qx,£lv) = fo(x,v), Vx,f E M 3 where SleG. 

The notion of asymptotic flatness is subtle. Intuitively one can think of an asymp¬ 
totically flat metric g a \, as a metric that far out tends to the Minkowski metric in some 
coordinate system. However this turns out to be a difficult definition to work with, 
because the coordinate invariance must be checked in all statements. It is not clear how 
to specify in which way the limits ”r —> oo” are to be taken in a coordinate independent 
manner. We will not give the precise definition of asymptotic flatness here, but just 
mention the fact that there are only a few symmetry classes compatible with asymptotic 
flatness where spherical symmetry and axial symmetry are the most interesting ones. 
The reason that cylindrically symmetric spacetimes are not compatible with asymptotic 
flatness is that the action group contains translation symmetry as mentioned in the in¬ 
troduction. Nevertheless, we can use another type of ” asymptotic flatness” condition 
because of the absence of translation symmetry in two directions of the 3-dimensional 
spacelike surface. Following ^SJ, we define asymptotically flat cylindrically symmetric 
initial data by: 

(i) (S, h^) is geodesically complete with S and S/G not compact. 
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(ii) There exists a G-invariant subset K of S such that h^ is flat on S\K. 

(iii) K/G is compact. 


(iv) /o = 0, in the complement of K. 


In j3j it is proved that the action of an abelian symmetry group together with this 
notion of asymptotic flatness and a reasonable energy condition, which Vlasov matter 
satisfies, exclude all possible spacelike surfaces and group actions except for S = M 3 with 
cylindrical symmetry or a periodic identification thereof along the translation axis. 

In this paper we are going to use coordinates t £ [0,oo), r £ [0,oo), z £ (— 00 , 00 ) 
and 6 £ [0, 2vr] which are adapted to the symmetry such that the two Killing vector fields 
are 


V 2 = 


A X 3 = ± 

dz 89 ' 

Hence the translation symmetry is along the z-axis and the rotational symmetry is 
around the z-axis. We consider a metric of the following form 


ds 2 = -ae^-^dt 2 + e 2 ^~ A dr 2 + e 27 (dz + Add) 2 + r 2 e~ 2 ^dd 2 , (4) 


where a, 17 , 7 and A depend on t and r. It should be remarked that this is not the 
most general asymptotically flat, cylindrically symmetric metric possible. However in 
the elctro-vacuum case it can be shown that a metric of this type is the most general. 
With this parametrization of the metric the Einstein-Vlasov system reads 
The Einstein-matter constraint equations: 


5=2™+_ ,/z^).,, 

r Zr z 


— = 7^ + — + 
r a 

Ot r 
2 r 


2 A 2 e 47 A 2 e 47 


4|->2 

= ae 2 ^~ A (Pi - p). 


+ 


4 or 2 


+ e 2 ^" 7 V, 


The Einstein-matter evolution equations: 


(5) 

( 6 ) 

(7) 


Ol-rr . ^tVt OL r . Q'r'Hr . 2 2 , 

m-Wrr = — + --- + —+ cn r - lt + 

+0 {A 2 - aA 2 ) + ae 2 ^) (P 2 - P 3 ), 


«t7t , «7 r , CXr'-fr i e 47 , 4 2 ^2\ , 

7tt - «7rr - 77W + ~ + “W” + 7772 (A ~ aA r) + 


+ 


2 a 
ae 2 ^ -7 ^ 


2 ' 2 r 2 

(/? - Pi + P 2 - P3) 


^ . A+at A r a r . . 

An — aA rr = — - 1 -—-4A*7t + 4aA r 7 r .-(- 

2 a 2 


aA r 


+ 2 are 2?? 47 S 23 . 


( 8 ) 


(9) 


( 10 ) 


4 















The Vlasov equation: 


9f + Jav 1 df 


dt 

+ 


dr 


2 ^)”° 


8 / 

dv 1 


+ 


(v*y 


r n 


- (m - 7 1) v 


i df_ 
dv 1 


a 7 r 


7t^ 


(u 3 )^ y/aAre 27 v 2 v 3 


+ 


df_ 

dv 1 


+ 


w 1 ?; 2 


s/aA r e 21 v l v 2 


+ 


r 

A t e 27 


+ V «-7r 


1 S 
V V 


df_ 

dv 2 

df_ 
dv 3 


+ 


(ID 


+ 


v 2 + 7 tv 3 


df 

-L- = 0 

dv 3 ’ 


where the variables v a are related to the canonical momentum variables p a by 
v° = \fae v ~ 1 p°, u 1 = e^-^p 1 , v 2 = e 7 p 2 + Ae 7 p 3 , v 3 = re~^p 3 , 
and the matter terms are defined by 

f(t,r,v)v°dv, 


P = 
J = 

Pk = 
S 23 = 


1 R 3 


IR 3 


f{t,r,v)v l dv, 


f (v k ) 2 

/ f(t,r,v)±-^-dv, k = 1,2,3, 

hi 3 v" 


f(t,r,v 


v 2 v 3 

1— h - 


dv. 


Jr 3 v- 

We assume that all particles have proper mass one, so 

9abP a p b = -1, 


and 


J° = \J 1 + (n 1 ) 2 + (v 2 ) 2 + (u 3 ) 2 . 


Observe that / does not depend on z and 6 because of the imposed symmetry as¬ 
sumptions. We assume initial data «o = a(0,r), 70 = 7 ( 0 , r), 71 = 7 ^( 0 , r), 770 = 77 (C), r), 
7 i = T 7 t( 0 , r), Aq = A( 0 ,r), A\ = A t (0,r) and /o = f(0,r,v) which is smooth and fur¬ 
thermore, according to the definition of asymptotic flatness we say that the initial data is 
asymptotically flat if there exists a set K which is compact in {f £ l 3 : z = const.}, such 
that /o = 0 in the complement of I\ and 70 , ? 7 o and Aq are constant in the complement 
of K. By a rescaling of z we can obtain 70 = 0 in the complement of K. It is important 
to point out that not only a plane is flat in two dimensions, also a cone is flat in two 
dimensions and the asymptotic behavior will in general be conical which admit ijq and 
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Ao to be different from zero in the complement of K. Furthermore 71 = 7/1 = A\ = 0 
in the complement of K. For the function ao we impose that lim a^fr) = C a , for 

r—>oc 

some constant C a . A natural choice is C a = 1. See |3 a for more details concerning the 
asymptotic flatness conditions. 

Regularity of the metric lead to boundary conditions of the metric on the axis of 
symmetry. In [3; the following lemma is proved. 

Lemma 1. Let A^ u be a smooth tensor field on the open ball B (e) C M 3 , e > 0, A^ IU 
invariant under rotation around the z-axis. Then there exist smooth functions x, '0; 

a, g and A invariant under rotations around the z-axis, such that A^ u takes the form 
(r,6,z — cylindrical coordinates) 

A flu dx tJ 'dx u = xdz 2 + ipr 2 dzdd + iprdzdr + ar 2 dr 2 + g (dr 2 + r 2 dd 2 ) + Xr 3 drdd. 

For a proof, see [U appendix C]. In our case this in particular implies 

A 0 (0) = 0(r 2 ),r -> 0, 
d o(0) = 0(1), r -»• 0. 

Furthermore the ratio between the circumference of a small circle around r = 0 and its 
radius must tend to 2tt as r —► 0. This and the fact that Ao(0) = 0(r 2 ) implies that 
770(0) = 0. Hence 

A 0 (0) =0(r 2 ),r—>0 , (12) 

do(0) = 0, r > 0. (13) 

Hence we call the initial data smooth, cylindrically symmetric and asymptotically flat if 
it is O 00 and satisfy the boundary conditions mentioned above, see equations d-CT 
with t = 0. 

From |3 p 90-91], we have that the symmetry assumptions on the initial surface are 
preserved by the time evolution. We will consider the solution given by Choquet-Bruhat 
0 which exists and is regular on some time interval [0, T). Since a is bounded on any 
closed subinterval and the characteristic system associated with the Vlasov equation, see 
(El and ®, yields that the matter is compactly supported on the same subinterval, 
it follows from m-m that the asymptotic flatness conditions hold for 7, 77, and A due 
to the hyperbolic nature of these equations. Hence the boundary conditions are 


lim f(t,r,v) = 0, 

r—>00 

(14) 

lim 7 (t, r) = 0, 

r—>00 

(15) 

lim ??(t,r) = C v , 

r—>00 

(16) 

lim Aft, r) = Ca , 

r—>00 

(17) 
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lim aft, r) = 1, 

r—>oo 

(18) 

A(t,r ) = 0(r 2 ),r —> 0, 

(19) 

o' 

T 

o' 

II 

"IT 

■+-T 

'l? 

(20) 


where the solution exists. To conclude this section we mention a lemma for cylindri- 
cally symmetric, asymptotically flat spacetimes, which show that these spacetimes are 
well-behaved in a certain sense. Namely, the Penrose singularity theorem m chapter 
9] says that the presence of a trapped surface implies the existence of a singularity, but 
the following lemma shows that trapped surfaces are absent. 

Lemma 2. Asymptotically flat, cylindrically symmetric initial data (S,g a p,K a p), with 
matter satisfying a reasonable energy condition, do not contain trapped surfaces which 
are either compact or invariant under the isometry group. 

For a proof and to see what is meant by a reasonable energy condition, see pj. 
The important thing for us is that Vlasov matter satisfies this condition. Thus possible 
singularities will not form due to the presence of trapped surfaces. 

3 The main theorems 

We are going to prove that if singularities form in the evolution, the first one occurs at 
the axis of symmetry. To obtain this result we will assume that the derivatives of the 
metric functions and the matter distribution function are well-behaved at the axis of 
symmetry so that the difficulties that arise there are eliminated. In the polarized case, 
i.e. when A = 0 however, we will not need any assumption on the metric due to the 
following theorem. 

Theorem 1. Given e > 0, let (R 3 , h^, K^, /o) be smooth, asymptotically flat, cylindri¬ 
cally symmetric initial data which satisfy the constraint equations and where the metric 
is coordinatized as ds 2 = —ae^^dt 2 + e 2 ^-^ dr 2 + e 2l dz 2 + r 2 e~ 2l d9 2 . Let ( f,g a b ) 
be the local smooth solution associated to the initial data and let [0, T) be the maximal 
time interval of existence. Assume that for r < e 

f(t,r,v) = 0,Vt £ [0,T),Vu £ R 3 . (21) 


Then T = oo. 

Due to the fact that the matter is supported away from the axis of symmetry one 
may suspect that we could use the vacuum results of 3] to control the field at the axis 
of symmetry. However this does not seem straightforward, due to the fact that a and 
at will be functions of t when r < e. Hence we can not reduce the nonpolarized Einstein 
equations to a wavemap from 2 + 1-dimensional Minkowski space to the hyperbolic plane 
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with the metric h tiu dx^dx v = d y 2 + du 2 in a vicinity of the axis of symmetry where 

the matter vanishes and then apply the methods of Christodoulou/Tahvildar-Zadeh. 
However as the theorem above shows we are able to remove the assumptions on the 
derivatives of the metric in the polarized case, i.e. A = 0. This is done by a combination 
of Sobolev estimates and the arguments in the proof of the following theorem. 

Theorem 2. Given e > 0, let (K 3 , h^, K^, /o) be smooth, asymptotically flat, cylindri- 
cally symmetric initial data which satisfy the constraint equations and where the metric 
is coordinatized as in @ . Let (/, gab) be the local smooth solution associated to the initial 
data and let [0, T) be the maximal time interval of existence. Assume that there exists a 
constant C e such that for r < e 


dh 1 


< C f , and 


d?A 


< C £ r, V/3 e N,Vt E [0, T), 


( 22 ) 


f(t,r,v) = 0,Vte[0,T). (23) 

Then T = oo. 

In assumption above we used the notation that f3 = is a multi-index 

and = d^dr 2 . 

Remark. By assuming that the support of the v 2 and v 3 momentum is uniformly bounded 
in time and that the modulus of “angular momentum” (which is more or less a product 
between the radial position of the particle and the velocity in ^-direction) and the velocity 
in the z-direction is bounded away from zero initially global existence is an immediate 
consequence of Theorem 1. This can be seen by using the Killing equations (EH) and cu 
together with the assumptions above to conclude that the matter never reaches r = 0 in 
finite time. 


4 Proof of theorem 2 

Given e > 0. If T = oo we are done so assume that T < oo. If we are able to obtain 
uniform bounds on all the functions (including their derivatives) we can extend the 
solution (/, g a b ) to t = T, and then apply the local existence theorem again with initial 
data at t = T, contradicting the fact that [0, T ) was the maximal time of existence. The 
outline of the proof follows [I] and we briefly summarize the strategy of the proof. In 
step 1 we bound the metric functions by using a conserved ’’energy”. In step 2 we bound 
Q{t) and the derivatives of the metric functions by using a light cone argument together 
with a crucial lemma, (Lemma 3). An important ingredient in the proof of this lemma is 
a rescaling of the momentum variables leading to cancellation of terms that are difficult 
to handle. In step 3 we bound the derivatives of the matter terms. Here the geodesic 
deviation equation plays an important role. It is then straightforward to see that the 
higher order derivatives can be bounded by similar arguments as will be outlined in step 
4. Whenever convenient we shall use the convention to write C for arbitrary constants 
even if they change from line to line. 
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Step 1 : (Bounds on the metric functions.) 

The asymptotic flatness condition m is crucial and leads to a uniform bound of r/ 
and a, because r/ r > 0. Indeed by © and m 


\v{t,r)\ = 


ij r (t , r')dr 


< 


poo 

/ rj r (t,r)dr 
Jo 


= lim \r](t,r)\ = C v . 


Hence |r/| is uniformly bounded on [0,T) x [0, oo). Note that the last equality above 
can be written 


poo pc 

/ rjr(t,r)dr = 

Jo Jo 


2 , It , A r e ^ , , ,,2(77—7) 


7r d-- + , 2 


+ 


4ar 2 


+ e 


dr = C V . (24) 


By the constraint equation ©. equation (HI and that P\ < p 


/ oo o poo 

— log a(t, r')dr' = — 2 J r'e 2 ^ -7 ^ (Pi — p) dr' < 

poo poo 

<2 / r' e Ari ~ A pd;r' < 2 / re 2 ^~ A pdr < 2C V . 

Jr Jo 

By the constraint equation © a r < 0. So lim a(t,r) = 1, implies that 

r—>00 

1 < a < e 2C ” =: (7. (25) 

Hence a is uniformly bounded on [0,T) x [0,00). 

By the asymptotic flatness condition m and the fact that |7(0,r)| = 0 in the 
complement of the set K, see m, we can choose R G (0,00), big enough such that 
|y(f,r)| = 0, V(f,r) G [0, T) x [R + Ct,oo). So by Holder’s inequality and equation (1271) 
we have for r > e 


|7(f,r)| = 


rR-\-Ct 

/ 7 r (t,r')dr' 

J r 


< 


rR+Ct 


2 / r R+Ct 


r'^dr' < C v 




R+Ct 


1 w 


-dr' < 
r 


< C v log 


R + Ct 


< C v log 


R + CT 


Hence 7 is uniformly bounded on [0,T) x [e,oo). If r < e then by Holder’s inequality 
and assumption Ha 


|7 (t,e) —y(t, r)| = 


^ 7 r (t,r')dr' 


< esup |7 r (t, -)| < C(e)e. 

r<e 
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Because 7 (t, e) is bounded due to the estimate above, 7 is uniformly bounded on 
[0, T) x [0, e). Hence 7 is uniformly bounded on [0,T) x [0,oo). 

Let R G (0,oo) be as above, then \A(t,r)\ = Ca , V(t,r) G [0, T) x [R + Ct, 00), 
which follows from the asymptotic flatness condition dm and that |H(0,r)| = Ca in the 
complement of the set K , see m- So by Holder’s inequality and equation EJ> 


I C A -A(t,r)\ = 


nR-\-Ct 

/ A r (t, r') dr' 
J r 


< 


rR+Ct 


< 


Ar'e 47 dr' 


rR+Ct a 2 4 j \ 2 


r -R+Ct 


<C V 


4 re 47 dr' < 


< C(e)(R + CT). 
Hence A is uniformly bounded on [0, T) X [0, 00 ). 


Step 2 : (Bounds on the first order derivatives of the metric and Q(t).) 
Define the support of the momentum Q(t) by 

Q(t ) := sup{|u| + 3 : 3(s,r) G [0 ,t] x [0,oo) such that f(s,r,v ) 7^ 0}, 

and define the two quadratic forms G and H by 

1 v 2 p 4 7 A 2 

G:= ^(— +7r) + ^(— + A 2 r ), 

2 a a 


(26) 




47 


a 4-^/ar 2 


-4/ . 


The aim is to show that the quantity 


T(f) := sup G(t, •) + Q 2 (t), 

r£[e,oo) 


is uniformly bounded on [0, T), by obtaining a Gronvall type of inequality 

r(t) < Cl + c 2 / r(s) log r( s )ds. 

Jto 

First of all, if <f> is a geodesic and X is a Killing vector held, then g(X, (j)') is conserved 
along the geodesic, see m p 442]. The particles follow the geodesics of spacetime with 
tangent p 7 so ( d z ) v and g^ w p^ ( do) u are conserved. Hence 
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9zzP z ( d z y + g ze p e (d z y = V 2 (t)e^’ R W = C, (27) 

9z6P z ( d e ) G + g ee p 6 ( d e f = V 2 (t)Ae^ R ^ + V 3 (t)R(t)e~^ R ^ = C. (28) 


Here R(t),V 2 (t) and V 3 (t) are solutions to the characteristic system associated to 
the Vlasov equation. The characteristic equation for R(t) is 


dR _ yJaV 1 
~ V° ’ 


(29) 


and we have especially that 


dR 


V^v 1 

ds 


v° 


< y/a < VC. 


(30) 


So R(t) is uniformly bounded on [0, T). More important though is that due to assumption 
d and that /o(r, v) = f(s, R(s, 0, r, v), V (s, 0, r, v)) it follows that 


R(t) > e > 0,vt e [o,r). 


(31) 


Furthermore the bounds on |y| and |H|, equation (TTT1) and that /o is compactly supported 
implies that |V 2 (f)| and |V 3 (t)| are uniformly bounded on [0, T) in view of equations 
d and d - So we can conclude that sup{|u 2 | + |u 3 | +3 : 3 (s, r, u 1 ) E [0, t] x [0, oo) x M 
with f(s,r,v ) y 0} is also uniformly bounded on [0, T). Hence in order to control Q(t) 
it is sufficient to control Q x (t) := snp{ |1 + 3:3 (s,r, v 2 ,v 3 ) E [0, t] x [0, oo) x R 2 with 
f(s, r, v) 7^ 0}. We also notice a few facts about the matter terms which follow from the 
discussion above. In the p — P\ term some cancellation occurs which is useful to take 
advantage of. Indeed we have 


..o K) 2 ' 


0 < (p - Pi) (t, r)= I ( v u - -f- ] f(t, r , v)dv = 

i + (v 2 ) 2 + (v 3 y" 


Ir 3 


lR3 \ V" y 

,P\ 2 +, 3 \ 2 \ 

f(t,r,v)dv < 


< / H ,( 1 + h 2 ) 2 ++) 2 ) i/i 


dv 


< C 


f 

J\v 1 \<Q 1 ( 


< 


1 + (u 1 ) 2 


(32) 


dv 1 


VI <QHt) ^l + (v 1 ) 2 
< C log Q\t). 


< 
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In a similar way the following estimates for Pf.(t,r) , k = 2,3 hold 


0 <P k 

< c 


(t,r) = [ 
Jh 


i?3 V 


(v k ) 

— 0 —f(t,r,v)dv < 


dv 1 


< 


|v 1 |<Q 1 (t) 1 + (v 1 ) 2 

< CTogQ 1 ^). 

The argument for S 23 is similar, hence 

S 23 < C log Q\t). 
Next, define two vector fields by 


( 33 ) 


(34) 


£ := - y/adr), 

C ■= dt + Vad r ). 

These vector helds are null vector fields, i.e. g a b^ a ^ b = 9 ab( a ( b = 0. Now by the 
evolution equations m and GDI 


t(G + H) = -j=(d t -yfaa r ) 


1.7 1 .9 e 47 At . . 2 

2 ( ^5 + 7 ’- > + &3 ( ^S + Ar) 


1 


3 

2a 2 


7t , W 

V2 k >/s 7r) 

e 47 2 . A t 

+ 8^72 ( 7^ + Ar 

8r 2 y /2 y/a \ 2a 

e 47 1 , 


7t«i . . 7t«r 

+ 7 «r — 7tr + -7^7 


V^T'Vr ^ 


+ 


A* 

'a 


^ + A tr -At r ) + 

2a 2 


+ i^T! ( 7l + Ar)2 ( 47t " 4v/ “ 7r + 

1 / It . y 
V / 2 ( \^ 7r) x 

e 47 1 /A 2 a r 2A r A t a r 2 A, 2 

+ -^— + * 
8^ 2 \/2 V at 


■rr J 

2 ^/a 


7ta r i a r 7 r f v / » 7 r- , V« e2(r? 7) A , 

1 ^ ,— 1 1 _ K ~r 


e 47 1 /A? 


+ 8r 2 i 


y.OL'f' 


+ 


2a ' 2\/a r 2 

_ + 4rA t e 2?? - 47 ,S 2 3 V 
V ar / 

+ 4v / ar-A r .e 2j? ~ 47 5 2 3^ = 


a 

2A t A r 


" r (G + H) + K£ 


>2(77-7) 


-C ( 7 ) + —C ( 7 ) + 


+ 


y/ 2 a 
S23C 2t] 
2 r 


Am 47 

C(A) + -^C(^), 

4- v /ar ,:> 


(35) 
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and 


Q(G-H) = -j={d t + V^d r ) 


-(^P= - 7r ) 2 + ^(A _ 

2 v \/a 8r 2 -i/a 


1 / It \ I IU 

= 7i ( 75~ 7r), ~ 


7« 7 tott , 7t«r 

T + 7tr - 7«r-X- 

a 2«2 2a 


V^OVr ^ "h 


p 4 7 9 4+ 

P 4 7 2 4 + 

_i_ c _ __c_2L _ 4 'i 

8r 2 V2 1 ^ rj 

e 47 1 , A t 


^L-^jL+A tr -A tr ) + 


« 2a2 

” ^A-r) + 

2-v/a 


+ “ A) 4 7t + 4^a 7 r - 


8r 2 -y/2 7/a 
1 / 7* 

\/2 V« 


/ 7t a r a r7r V«7r . \/ae 2(,) 7) \ 

VS + ~ + -2- K J + 

+ |L-L (-47 + _ M + 4,4,7" <"S a ) + 

8r A V2 V a! ct Var / 


+ 


7 


A 2 a r 2A t A r 


s 47 1 


8r 2 ypi \ -y/a 

{G-H) + Ke 


- 4\/arA r e 2T7 47 = 


2 ( 77 — 7 ) 


-f (7) + — f (7) + 


+ 


\/2a 

S 23 e 2ri 


Ha)-^ALha), 


(36) 


2r 71 v ' 4-y/ar 3 

where k = p — P\ + P 2 — P 3 . Let f±(s) be the unique solutions to the equations 


AA = ±v/a(s,f ± (s)), f±(t) = r , 

which exist due to the theory of ODE because a is a smooth function and a > 1 by (ESI). 
Note that because r > e we can find a largest time to, 


to < T 



(37) 


such that if we go backwards along the integral curve for £ we do not come closer to 
the axis than r = Observe also that we can control the different functions on the 
hypersurface t = 1 0 due to the local existence theorem of Choquet-Bruhat [3 • Integrate 
the equations (1351) and El for the quadratic forms above along the integral curves 
corresponding to the vector fields ( and £, from the hypersurface t = to < T, to the 
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point ( tr ) G {to,T) x [e,oo) to obtain 


(G + H) (t,r) — (G + H) (to,f-(t 0 )) = 


rt / ~ 2(^-7) S^e 2 ^ 


Olf 


a 


(G + tf) + 


Ace 


-C (7) + 


2 r 


C(^) (s,r_(s))ds+ 


+ 




+ 7 - 7 = 3 C (-4) ) (8,f-(s))ds, 
4 ^/ar 5 


and 


(G-if) (t,r)-(G-H) (t 0 ,M*o)) = 


^ (G — H) + Ke 


2('<?— 7 ) 


+ 


■'to 

f 

1 to 


a 


< (7) + 


5 2 3e 2r? 


2 r 


£ (-4) (s,f+(s))ds+ 


—C (7) - 7-7=^^ (^ 4 ) ) (s, r+(s)) ds. 

r 4yar J 


Add the equations to obtain 


G(t,r ) = 

= 7 [G(i 0 ,M*o)) + G(t 0 ,f + (t 0 ))] + 


+ 7 [#(*(), M*o)) - H(t 0 ,r + (t 0 ))] + 


+ 


+ 


+ 


+ ^> 


*0 


'to 


'to 


to 


, r a (G + H) + —j— C ( 7 ) + -^_C (^4) (A r_ (s)) ds+ 


— C (7) + 7T7=-^C (- 4 )) ( 8 ,f-(s))ds+ 

r 4 ycur / 

(G - fl) + ( 7 ) + (A)] (s, f + (<0) *+ 


a 


—£ (7) - 3 C (^4) ) (a, Ms)) ds - 

r 4\/a;r' :> 


Observe that G + H is always nonnegative and also from the constraint equation © 
that < 0 , so the term -^= (G + H) can be dropped which simplifies the estimate 
but it is not necessary. Observe also that H < G because of the elementary inequality 
2 ab < o 2 + b 2 . Take the supremum with respect to r 6 [e, 00 ) and use that the matter 
is compactly supported, equations (EH)-imp, together with the constraint equation 0 , 
assumption EH and observe that the radial distance from the axis is bounded from 
below by then 
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sup G(t,-)< sup G(t 0 ,-)+ sup H(t 0 ,-) + 

rG[e,oo) rG[e,oo) rE[e,oo) 

/■* ( Ke 2 ^ S 23 e 2 A rp^r - -r , 

+ SU P / G ---1--— \/G(s, ■)d-s+ 

r£[l,oo)Jt 0 \ ^ ^ / 

+ sup f Cr\foi^ r] ~'^ log Q 1 (s) (G(s, •) — H(s, -))ds+ 
r£[ |,oo) Jto 


+ sup 


r£[f,oo) z Jt 0 


/’(^ (7) + 


4^/ar 3 


C (^4) ) ( s )') ds+ 


+ sup 


1 


rt 


rG[§,oo) ^ Jto 


( 7 ) - 77 =^ (^)) 0> •) ds < 

r 4 y/otr* 


(38) 


<Ci+C 2 (e) f log Q 1 (s) sup G(s,-)ds + C 3 (e) f sup G(s,-)ds. 

Jto rE[e,oo) Jto rG[e,oo) 

Because we are interested in the case when sup G is large, we assume throughout the 

rG[e,oo) 

paper that sup \/G < sup G. The following crucial lemma will yield a bound on 

rE[e,oo) rG[e,oo) 

|Q 1 (t)| 2 in terms of G. 

Lemma 3. Let Q 1 and G, be defined as above, then 


{Q\t )) 2 <c 1 + c 2 


'to 


(Q 1 (s))~ logQ 1 (s) + sup G(s, 

r*G[c, 00) 


ds. 


Proof. Define 17 1 (s) := e 71 7 D 1 (s). 

The characteristic equation for D 1 (s) reads 


(Or - 7-) + v° - (v, - 7.) r'+ 


+ 


(o: 


« 7r 


1\ iy> sJaA r e 21 V 2 V 3 


R V° 


+ 


R D° 
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Hence 


^-{U\s)f = 2U ldUl 


ds 


ds 

= 2U 1 e 71 ~ 1 
= 2U 1 e 11 ~ 1 
+ 2U 1 e v ~ 1 

+ 2U 1 e v ~' y 


, , T . , dR dV 1 

(lit ~7t)V +(Vr-7r)V ^ + ^ 


{lit - ItW 1 + ( T) r - J r )V 


rifaV 1 


+ 


U° 

^ V° —fr* — 


- [V® {Vr - Ir) + v° - (rjt - 7 1 ) V 1 


2 ^ J 


+ 


2\ 2 


U° 


a 7 r 


1 \ (U 3 ) 


3t 2 


R U° 


+ 


I 2U 3 c n ~ ^ Are2lV2V3 
+ ~ U R U° 


= 2 U 


i' 


+ 2 U 1 - 


V° 

0 * 1-1 

“u^ 


- ( \fa {p r — 7 r ) + (! + (^ 2 ) 2 + {V 3 Y) 


+ 


^7,(^) 2 -v^(7,-^)(F 3 ) 


3\2 


+ 


ie v 7 / s/aA r e 21 1 


+ 2U V° \ R 
= T 1 + T 2 + T 3 , 


V 2 V 3 - 


OLj' 


C v l f 


2 \[ol 


where 


T\ = 2 U 1 


T 2 = 2 U 1 - 


0 * 1-1 


u° 


0 * 1-1 


u° 


- f Va{dr -7 r) + 


Or \ 

2 ^EJ 


(1 + (V 2 ) 2 + (U 3 ) 2 ) 


^7r(^ 2 ) (^ 3 ) + 


2 \faA r e 


27 


R 


-U 2 U 3 


t 3 = -c/ 1 - 


0 * 1-1 


Otr * 


1\2 


U° ^ V ) ' 

Use the constraint equations © and © then T\ becomes 


\Ti\ = 


2 U 1 


2 U 1 


0 * 1-1 


U° 


0 * 1-1 


-[Va(Vr- lr) + 


Cfc r 


2 y/a. 


(l + (U 2 ) 2 + (U 3 ) 2 ) 


u° 

0*1-1 r 


— y/oL I i?7r 4-—I" 


R'y 2 H 2 e 47 A 2 e Al \ 


a 


4 R 


+ 


4 Ra 


■ 


(1 + (U 2 ) 2 + (U 3 ) 2 ) + 


+ 2 U 1 -^- [-V^Re^^p + ^7r - RV^e 2(v ~ l] {Pi ~ p) (l + (V^ 2 ) 2 + (U 3 ) 2 ) 
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2 U 1 


0^—7 


yo 

=57-7 


2 , ^7? , 7 l 2 e 47 , ^?e 47 


— \/Ol I R'Jr H-h 


a 


47? ARa 


+^--^ui+(vr+(v 3 ) 2 ) 


+ 


- 2t7 1 —Q-i?v / ae 2 ^ ?_7 - ) Pi (l + (H 2 ) 2 + ( V 3 ) 2 ) 


y° 

Notice that from step 1 we have bounds on |r/| , |7| , a and |A| and recall that |T^ 2 (i)| 
and |l/ 3 (f)| are controlled and e < R(t) < C. The matter terms can be bounded 
terms of Q 1 (t) because 


in 


Pi 


(t,r) < p(t,r) = [ f(t,r,v)v°dv<C\\f 0 \\[ [v 1 \ dv 1 < C (Q 1 (t )) 2 

Jr 3 ■'b 1 |<Q 1 


Furthermore notice that 


XP_ 

yo 


< C. 


(39) 


(40) 


Hence 


|Ti| < Ci sup G(s,-) + C 2 (Q 1 (s))' 

rS[e,oo) 


For T 2 use CD to obtain 


|T 2 | = 2 U 1 ^- |^ 7r (V^ 2 ) 2 - VS ( 7r - ^ ) (y 3 Y + 
<C sup G(s,-). 

re[e, 00) 


i\ (V 3y + ^ 2 \* 

K / Jr i 


H z l/ 3 


< 


By using the constraint equation 0 and equation dS2D 


|T 3 |= 2U l ' 


=»?-7 




H° 2^a 




1\2 


=»?-7 


2C/ 1 — r R^e 2 ^ (P 1 - p) (H 1 ) 2 


V° 


< 


< fiU'e^Ryfae 2 ^ logQ 1 ^)) I!/ 1 ] < C (Q 4 (s)) 2 logQ 1 ^). 
Add the estimates for |Ti|, |T 2 | , and |T 3 | to obtain 
d 


ds 


(^(•s)) 2 < Ci sup G(s,-) + C 2 {Q 1 {s)Y log Q 1 (s). 


rS[e,oo) 


Integrate the inequality above, note that H 1 (t) = e ^ 7 )£/' 1 (t) and that \r]\ and |^y|, 
are bounded. Then take the supremum of |H 1 (t)| to obtain 


(Q 1 ^)) 2 <Ci + C 2 / 

Jt ( 

and the lemma is proved. 


to 


(Q 1 {s )) 2 logQ 1 (s) + sup G(s, •) 

rE[e,oo) 


ds, 


□ 
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Remark. It is necessary to rescale the momentum V 1 to obtain cancellation of derivatives 
of the metric in the characteristic equation for V 1 . In [T] this was not necessary. The 
important cancellation in [T] was instead obtained by a careful analysis of the character¬ 
istic equation by observing the fact that the terms had the “right” sign. However due to 
the fact that we have control over the metric functions 7 and 7 this rescaled momentum 
is essentially the same as V 1 . 

By the estimate HI for sup G(s,-) and Lemma 3 

rg[e,oo) 

r(t) <Ci + c 2 f r(s)io g r(s)ds,te[t 0 ,T), 

Jt 0 


which is a Gronvall type of inequality. Because of the local existence theorem we have 
control over all quantities in the interval [0, to]- Thus sup G(s,-) and Q 1 (t) are uniformly 

rE[e, 00) 

bounded on [0, T). Hence it follows that the matter terms are also uniformly bounded 
on [0, T) in view of equation (ETUI) above. Note that due to assumption C2) we have 

control over sup G(s,-), 
re[ 0 ,e) 


and 


l,7f 2 \ e4 \^t ,2 

~ — +7r ! ) + 5"2 — +A V 
2 a 8 r- a 


G< sup G(t, •), 

re[ 0 ,oo) 


where every term on the left-hand side is nonnegative. This estimate and that a < C by 
P51) . implies that |7t| and |7 r | are bounded when r e [0,oo) and that \A t \ and \A r \ are 
bounded on compact intervals [0,1?]. Recall that by the asymptotic flatness conditions, 
\At\ and \A r \ vanishes when r > R + Ct , with R and G as before, see the discussion after 
equation (HJ, so that \At\ and \A r \ are bounded everywhere. 

The bound on |a r | follows from the constraint equation © and that a r = 0 when 
r > R + Ct, by the asymptotic flatness condition. Hence 


Gif 


2 rae 2 ^ ^ (P\ — p) 


< GlogQ 1 ^). 


Prom the constraint equations © and © the estimates on |r/t| and rj r follows. Fur¬ 
thermore | a* | < ClogQ 1 (t) < C, see Lemma 6 in the appendix. 

Step 3 : (Bounds on the first order derivatives of the matter terms.) 

To obtain bounds on the first order derivatives of /, it is sufficient to obtain bounds 
on dR and dV because the solution / can be written in the form, 


f(t, r, v) = /o(R(0, t, r, v), V (0, t , r, v)), 

where ( R(s , t, r,v),V ( s, t, r, v )) is the solution to the characteristic system associated to 
the Vlasov equation, see m, with R(t,t,r,v ) = r and V(t,t,r,v) = v. Note that a 
partial derivative for / reads 


dj_ = dfodR dfodV 
dz dr dz dv dz ’ 
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where z denotes either t, r or v. This part of the proof is almost identical to the proof 
of Lemma 3 in [T] but we have to use assumption CS1) that the matter never reaches the 
axis in finite time. 

Lemma 4. Let R(s) = R(s,t,r,v ) and V k (s ) = V k (s,t,r,v), k = 1,2,3 be a solution 
to the characteristic system associated to the Vlasov equation. Let d denote dt,d r or d v 
and define 


T = dR, 


Z 1 



7tV° (V 0 ) 2 - (V 1 ) 2 + (V 2 )' 2 - (V 3 ) 2 
Va ( V °) 2 - (F 1 ) 2 


/ A r e 21 LV 2 L 3 _ A t e 2 ' 7 V°V 2 V 3 \ 

V R (F°) 2 -(F 1 ) 2 V^R (F 0 ) 2 - (V 1 ) 2 ) + 


V 1 f(v 2 ) 2 - (F 3 ) 2 ) 

7 r -- -5-5 — -dR 

(V°) - (F 1 ) 2 


dR+ 


Z 2 = 8V 2 + 7 r V 2 dR, 

^ dR. 

Then there is a matrix A = (ai m ) ,l,m = 0,1, 2,3 such that 

(T, Z 1 , Z 2 , Z 3 ) T , 

satisfies 

f = (41) 

and the matrix elements are all uniformly bounded on [0, T). 

The proof is obtained by a long calculation, see [TJ for a sketch of the proof. 

By standard results from ordinary differential equations Q is uniformly bounded on 
[0, T), and hence the derivatives of the matter terms are uniformly bounded on [0, T). 
It is then immediate to obtain bounds on dR and dV l ,i = 1, 2, 3. 

Step 4 : (Bounds on second order and higher order derivatives.) 

This section also follows T and the methods used in step 2. Hence derivatives of 
arbitrary order both for the field and the matter can be uniformly bounded on [0, T), so 
the solution can be extended beyond T which is a contradiction. Hence T = oo and the 
proof is complete. 


Z 3 = 8V 3 + 


A r e 2 W 2 

~R 


-7rF d 
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5 Proof of theorem 1 

Proof. Given e > 0. If T = oo we are done so assume that T < oo. If we are able to 
obtain uniform bounds on all the functions (including their derivatives) we can extend 
the solution ( f,g a b ) to t = T, and then apply the local existence theorem again with 
initial data at t = T, contradicting the fact that [0, T) was the maximal time of existence. 
To exploit assumption m we divide ([0,T) x M. 3 ,g a b) into an interior-, r < |, and an 
exterior region, where r > |. The bounds in the exterior region rely on the proof of 
Theorem 2. However in that proof the assumption (E3 that the derivatives of 7 are 
controlled in the interior region was used in a crucial way. By an examination of the 
light-cone argument in the proof of Theorem 2 we see that it is necessary to obtain an 
estimate of the type (cbyl < Q l (t) in the interior region, and then apply the light cone 
argument as in the proof of Theorem 2 to the exterior region. To obtain this “interior” 
estimate we use a standard Sobolev method. 

From the proof of Theorem 2 we have that a and |?y| are uniformly bounded on 
[0, T) x [0,oo). Furthermore |7| is uniformly bounded on [0, T) X [e,oo), for any e > 0. 
From m and C9> we have 


V 2 {t)e l[tMt)) = C , 
V 3 (t)R{t)e~^ tMt)) = C. 


The discussion before and after equation m yields that 

\v 2 (t)\<c, 

\V 3 (t)\ < c. 


Hence 


sup{|u 2 | + |u 3 | + 3 : 3(s, r, u 1 ) G [0, t) x [0, oo) x (—oo, oo) with f(s, r, v) / 0}, (42) 

is also uniformly bounded on [0, T). So in order to control Q(t), see equation (l26l) . it is 
enough to control Q 3 (t) := sup{|u 1 |+3 : 3(s,r,v 2 ,v 3 ) G [0, f)x[0,oo)xR 2 with f(s,r,v) / 
0}. By the constraint equation 0 and assumption m 

a r {t , r) = 0, t G [0, T),r < e. (43) 

Hence a = aft, r ) is independent of r when r < e. Observe from the discussion above 
that assumption en> implies that the matter terms vanish when r < e. 

Now 7 is a rotationally symmetric solution of 

Dn'y = F, (44) 


where 


□tf = - 


d± 

Oi 


OLt 


CXnr 


(Xq, 


+ di + d, + -^dt + ~^d x + 


2a 2 


2a 


2a 


J yi 


(45) 
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with x = r cos 9, y = r sin 6 , and 


P^iv-'y) 

F = - —( p -P 1 + P 2 -P 3 ). (46) 


Let 



(47) 


Define the backwards truncated “cone” 0 = {(t,x,y) : to < t < T, \J x 2 + y 2 < f(t)}, 
where f (s) is the unique solution to 


rj p _ 

— = —yja(s, r(s)), f{T) = 5. 

Define 


(48) 


B t :={{x,y) : y/x 2 + y 2 <T + 5 - t}, (49) 

and let Ct 0} t be the surface of the cone. Choose to big enough such that the radius 
of Bt 0 is less than 25. This is always possible, and in view of the local existence theorem 
due to Choquet-Bruhat (Sj we can control all quantities on the hypersurface t = to- Then 
we can use the following special case of lemma 2 from El- 

Lemma 5. Let u G Hf oc (M . n ) with Bs(x) = {y G M n : |y — x\ < 5} on arbitrary ball. 
Then with |x|, 

2 

\u{x)\ < c^5 k ~ l || 3 k u \\ l2 {§s{6)) . (50) 

k =0 

Remark. Define fin ■= \foidxdy as the measure associated to the metric Hij. That 
\Ja > 1 in our case is immediate from the constraint equation 0 , see (ESI). 

Hence 

2 

\dl{t,x,y)\ <^5^ II dV dl (Vr) || L 2 . (51) 

u =0 

A typical term to be estimated in the Sobolev inequality above reads 

II d u d'y (i,y) ||(52) 

for 0 < v < 2. Observe that d u = d^dy 2 where iq +17 = v. Note also that ^ (p 2 + 07^ + 
a Ty) ~ yJ&jTtTx ~ > 0 and apply the local energy identity which reads 

I! Wl2(B t , m )= (53) 

= 11 d v d-i Hia (B - [ [^( 7 1 + «7x + «7 y) - Va-ltlx - \foi-Ttl y \da+ (54) 

•ICtn.t z r r 
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— [ + a lx + a'j%)y/adxdydT — [ d u Fd Ur jtV^dxdydT < (55) 

2 Jn 2 « Jn 

<11 d v d 1 \\l 2 { B t u H ) + c sup |a t (t,v)l / II ^7 II!2 ( b t ,^) dr, (56) 

(x,y)eK 2 7t 0 

where da is the induced measure on the surface of the cone. 

So we obtain a Gronvall inequality for || d u d'y |^ 2( ' Bf jlH y Hence 

(t—to) sup |o!i 

II ff'djit, -, ■) ||| 2(St ^)< c ie (57) 

From Lemma 6 in the appendix we get that \at\ can be estimated by 

\a t (t,x,y)\ < C 2 log Q 1 (t). (58) 

Finally 

II 9"07(V, •) lli. (a , M) < c, [e 1 («)] 02(, - , »>. (59) 

From Lemma 6 in appendix we see that C -2 do not depend on to, so by choosing to such 
that C*2(t — to) < 1 and of course such that the condition that the radius of Bt 0 is less 
than 25 holds, then 


\dl{t,x,y)\ 2 < CQ 1 {t),t G [t 0 ,T), ^ x 2 + y 2 < 8 . (60) 

Now we can rely on the proof of Theorem 2 with A = 0 and with e replaced by 
| in that proof. Furthermore observe that we have to replace assumption © by the 
estimate 

sup G(t,-) < CQ l (t), (61) 

re [f’f] 

which was obtained above. Recall the definition of G(t, r ) and H ( t , r ) 


G:=H^+7. 2 ). 

(62) 

2 V« J 

\/a 

(63) 


By using ini instead of assumption ©, after choosing to as in the discussion after 
equation dm and of course such that the radius of B to is less than 26 holds and also 
that m hold. Then the light cone estimate in the proof of Theorem 2 becomes, see 
equation dm- 


sup G (t, •) < 

re[|,oo) 

sup G(to,-)+ sup H(to,-) + 

re[f,°o) re[f,oo) 
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+ sup / 

re[|,oo) Jto 


c 


Ke 


■\/G{s, -)ds+ sup / Cry/ae 2 ^ 1 log Q 1 (s) (G(s,-) — H(s,-)) ds+ 


re[|,oo) Jto 


+ 


c 

sup / \fa — G(s, -)ds 

:|L nrd Jto T 


re [f’°°) 


< Ci + C 2 (e) / logQ^s) sup G(s, -)ds + C 3 (e) / sup G(s,-)ds+ 


'to 


re[|,oo) 


'to r&[j ,oo) 


C 4 (e) I logQ 1 (s) sup G(s,-)ds + Cs(e) f sup G(s,-)ds< 
Jto re [|>l) Jto re[|,|) 


<Ci + C 2 (e) f log Q l (s) sup G(s, -)ds + C 3 (e) f sup G(s,-)ds+ 

Jt 0 re[|,oo) Jto re[f,oo) 

+C 4 (e) f Q 1 (s) log Q 1 (s)ds + C 5 (e) f Q 1 (s)ds. 

■Jt 0 Jt 0 

By combining this with Lemma 3, which holds without any changes, a Gronvall type 
inequality is obtained for 

,2 


It reads 


T(s) := sup G(t,-) + (Q 1 ^))' 

re[f ,oo) 


r(t) < Ci + c 2 [ r( s ) log r( s )ds. 

Jto 


(64) 


(65) 


Hence T(t) is uniformly bounded on [to, T) and as mentioned above we have control over 
all quantities when t £ [0, to]- So T(f) is uniformly bounded on [0, T). Together with 
equation CT this implies that G is uniformly bounded on [0, T) x [0,oo). 

To conclude we obtain uniform bounds on the metric functions and all the first order 
derivatives and on the matter terms, see the proof of Theorem 2 and the appendix. Now 
we can continue in the same way as in the proof of Theorem 2 with the only modification 
that e is replaced by | and that assumption (PI) is replaced by 


3^7 <CQ\t)yfJ> l,te [t 0 ,T),^/x 2 + y 2 <6, 


( 66 ) 


where Q l (t) < C in view of the estimates above because a is independent of r when 
r < e. Hence Theorem 2 is proved. 

□ 

Remark. In the nonpolarized case this Sobolev method does not work. This is to be ex¬ 
pected in view of the deep proof of Christodoulou and Tahvildar-Zadeh on the spherically 
symmetric wave maps from flat space into a class of Riemannian manifolds. 
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6 Appendix 

Lemma 6. With all the functions as above the following estimate holds 

\a t (t,r) | < ClogQ 1 ^),* € [0,T), (67) 

where the constant C only depend on the maximal time of existence T. 

Proof. Observe that we have uniform bounds on a and \p\ when (t,r) £ [0 ,T) x [0, oo) 
and a uniform bound on |7| when (t, r ) £ [0, T) x [e, oo), see step 1 in the proof of Theorem 
2. Furthermore assumption m implies that |U 2 (t)| and |U 3 (i)| are uniformly bounded 
in time, see step 2 in the proof of Theorem 2. Now integrate the constraint equation © 

at j t,r ^ = 2 re 2 fo-7)(p p ). 

a(t,r) V 

Use the asymptotic flatness condition m 

lim a(t, r) = 1, 


to get 


Hence 


/ OO 

r , e 2 ( ?? - 7 )(p — P\) dr'. 


aft , r) 


= 2 


r'e 2(r? 7) [2 (rj t - 7 t )(p ~ Pi) + d t (p - Pi)] dr' = 


= 2 


POO 

/ re 2(r?_7) [2 (rjt - 7t)(p - Pi) + <9 t (p - Pi)] dr. 

Jo 


Observe that assumption m implies that at = at (t), r < e. The first term in the 
integral above can be estimated as follows 


4 / re 2(??_7) (r] t - 7 1 ) (p - Pi)dr 


<C\ogQ l (t) j r\r] t - 7*1 dr < 


<C\ogQ l {t) j r (|r?i| + |7* |) dr. 

Jo 

From the constraint equations m and © for 7 

\Vt\ < CV 

Recall that the derivatives of the metric have compact support so by Holder’s inequality 

POO POO 

C log Q l (t) l r (\r) t \ + \^t\) dr < ClogQ 1 ^) / r fq r + |7 t |) dr 
Jo Jo 


< 
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< Cl log Q 


r 

Jo 


< ClogQ 1 (t) 


poo 

ilrdr + C 2 \ogQ 1 (t) / r|7, 

Jo 

1 / 

poo 

/ l7i 
J o 


'tl dr < 


1 1 2 rdr 


r/ r dr + 

< ClogQ 1 (t), 


in 

R+Ct \ 2 


rdr 


< 


where we used that 


poo poo 

/ I'Yt | 2 r-dr* < / r\ r dr < C. 

JO JO 


Hence 


poo 

4 / re 2(r?_7) (rjt - it) (p - Pi)dr 

Jo 


< ClogQ 1 ^) 


For the second term in the integral expression of at above use the Vlasov equation to 
obtain 


poo poo p 

\ / re 2 {jl ~^dt{p — P\)dr = 2 / 

Jo Jo J i 


re 


2(^—7) 


■) 1 + O 2 ) 2 + (r 3 ) 2 r 


+ (^/a(p r - 7 r) + 


Gif 


0 JR 3 


id 9 / 


+ 7tu + v ct 


2y / a / Sr 1 

rVh df 


(v 2 ) 2 !— 1 (r 3 ) 2 1 \ 9/ 

^r— + Mir ~ ~) — + (7t - 7t)« J ^1 + 


0 I Qy2 


—2 + V^(-- 7 r) 


1 S 
V V 


df 


dvdr 


We estimate different terms separately as follows. By the conserved quantity 01, the 
observation in the beginning of appendix and partial integration 


r ,■ 

Jo J R : 


2^)1 + (r 2 ) 2 + (r 3 ) 2 Mdf, , 
re w v -n- \Ja—rr—dvdr 

7U y\J 


yV 


2 re 2(v ~ 7) 


m[ 

J 1 : 


1 + (r 2 ) 2 + (v 3 ) 2 V 1 


: f(t,r,v)dv+ 


+2 


poo p 

JO J R ; 


1 + {v 2 ) 2 + {v 3 ) 2 v 1 


J(t,r,v ) 


e 2(r? - 7) V« 


+r[e 2( ' v ^M^ilr ~ lr) + re 2(v ^ M=\ 

2y a 


dvdr 


< 


< Ci log Q 


poo p 

\t) + 2 / 

Jo Jr- 


1 + (r 2 ) 2 + (r 3 ) 2 v 1 


poo 

CXf' 

to 

2 v / a 


C 2 log Q 


rc 

\t) 

Jo 


re 


2(^—7) 


■f(t,r,v)r'e 2( ' ri 7 -*— '-=dvdr' < 

2 V« 


dr < Ci logQ 1 (^) + 


pdr < 
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Clog Q 1 ^). 


By partial integration 


re 


2(77—7) 


1 + (u 2 ) 2 + (v 3 ) 2 


Jo J R 3 
Again by partial integration 




\ 0 df 

= V — - T 


- 7r) + x-^= v^—^jdvdr 
2^/a J ov L 


= 0. 


/ 0 J R 3 


/ 

it 


re 


2(77—7) 


1 + (u 2 ) 2 + (u 3 ) 2 


2\2 


fa 2 ) 


1 (r 3 ) 2 \ df 


-Vair— n- + \/«(7r - -) 


r u° / Sr 1 


-du<ir+ 


+2 


2 ( 77 — 7) 1 + fa ) + fa 3 ) 7 N 1 df 

re 1 77 --(■ r] t - 7t) u -^dvdr 


' J q f R3 re2(V 7) /(^ r ^) (! + (v 2 ) 2 + (v 3 ) 2 ) (-\/« 7 r (^ 2 ) 2 ^4 + 


1 277 1 / 1 ( 7 , 1)2 \ 

+ A(7r - -)(^ 3 ) 2 -^0)4 + - 7t) + ^0)3 ) dV(ir 


< C log Q\t). 


The last terms can be estimated in a similar way by partial integration. Sum all estimates 
to obtain 

\&t(t,r)\ < ClogQ 1 ^). 

Furthermore it is immediate consequence from the proof above that C only depends on 
the maximal time of existence T. □ 


Acknowledgement: I want to thank Hakan Andreasson for suggesting the problem 
and for useful discussions and comments. 
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